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TECHNICAL NOTE 2267 

INELASTIC COLUMN BEHAVIOR 
By John E. Duberg and Thomas W. Wilder, III 


SUMMARY 


The significant findings of a theoretical study of column behavior 
in the plastic stress range are presented. When the behavior of a 
straight column is regarded as the limiting behavior of an imperfect 
column as the initial imperfection (lack of straightness) approaches 
zero, the departure from the straight configuration occurs at the 
tangent -modulus load. Without such a concept of the behavior of a 
straight column, one is led to the unrealistic conclusion that lateral 
deflection of the column can begin at any load between the tangent- 
modulus value and the Euler load, based on the original elastic modulus . 

The behavior of a column with vanishing initial lack of straightness 
at loads beyond the tangent-modulus load depends upon the stress-strain 
curve for the material. A family of curves showing load against lateral 
deflection is presented for idealized H-section columns of various 
lengths and of various materials that have a systematic variation of 
their stress-strain curves. These curve s show that, for columns in which 
the material stress-strain curves depart gradually from the initial 
elastic slope as is characteristic of stainless steels, the maximum column 
loads may be significantly above the tangent -modulus load. If the depar- 
ture from the elastic curve is more abrupt, such as for the high-strength 
aluminum or magnesium alloys, the maximum load is only slightly above the 
tangent -modulus load. 


INTRODUCTION 


Until recently, the double -modulus theory had generally been 
accepted as the correct theory of column failure in the inelastic range 
of stress. This 'theory, originally developed by Considere and Engesser 
and later extended by Von Karman (see reference 1 for a discussion of 
this development), predicts that the load at which bending starts and 
the maximum load that a pin-ended column can support are the same and can 
be obtained from the Euler equation 


PE = 


j* 2 EI 
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by substitution of a reduced modulus for Young's modulus. (The symbols 
used in this paper are presented in appendix A. ) The reduced modulus 
is obtained by assuming that, at the start of bending of an originally 
straight column, the direction of straining of the elements on the con- 
vex side of the column reverses. When such reversal of strain occurs 
in the plastic range, increments of stress are related to increments of 
strain by the original elastic modulus i The part of the cross section 
over -which the strains reverse is determined by the condition that there 
shall be no change in load during the bending process. 

In I 9 V 7 , Shanley (reference 2) was able to show for a simplified 
column that, if the load is allowed to increase during bending, bending 
of the column can start at a lower load than the reduced -modulus load. 
The load for which he showed this to be true was the tangent -modulus 
load, which can be obtained from the Euler equation by substitution of 
the tangent modulus in place of Young's modulus. Shanley drew conclu- 
sions concerning the behavior of columns on the basis of the behavior 
of the simple model and certain experimental observations. 

In order to clarify the behavior of columns in the plastic range, 
a theoretical study was made with the following threefold purpose: 

(1) To establish the load at which a column starts to deflect, 
designated the critical load in this paper 

(2) To study the mechanism of column action beyond the critical 

load 

(3) To establish the relation between the maximum column load and 
the stress-strain curve for the material 

In order to make this study, two models were chosen: One model, similar 

to Shanley' s, was a spring-supported rigid column - that Is, one that 
had a concentrated flexibility - and the other, an idealized H-section 
column that had its flexibility distributed along its length and con- 
sisted of two concentrated flanges separated by a web of negligible area 
but of infinite shear rigidity. 

The significant results of this study were presented without proof 
at the Structures Session, Eighteenth Annual Meeting of the Institute of 
the Aeronautical Sciences held January 23-26, 1950, in New York, and 
were subsequently published by that organization as reference 3 . The 
purpose of the present paper is to give the details of the analysis 
which were not included in reference 3. Since the original presenta- 
tion, several other investigators have published the results of their 
researches on inelastic column action. (See references ^ to 6 . ) The 
results in these references are in substantial agreement with those 
obtained in this paper. 
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THE CRITICAL LOAD 


In the elastic range of stress the critical load for a straight 
column - which must again be emphasized as having the restricted meaning 
in this paper as the load at which bending starts - is unique and is 
given by the well-known Euler formula. The situation is not so simple 
in the inelastic range of stress, however, and the source of the diffi- 
culty lies in the character of the stress-strain relations in the plastic 
range. In the Inelastic range of stress, at least for uniaxial states 
of stress, increments in stress are related to increments in strain by 
the tangent modulus of the material, but decreases in stress are 
related to strain by the original elastic modulus. 

In order to demonstrate the lack of uniqueness in defining a load 
at which bending starts in the inelastic range, the simple spring- 
supported column shown in figure 1 was analyzed. Details of the analysis 
are presented in appendix B. The inverted tee is rigid and free to 
rotate and translate vertically about the intersection of the cross of 
the tee. The column is supported at each end of the horizontal leg by 
identical sets of two elastic springs . One spring is fixed at the far 
end; the other spring has a ratchet attached to its far end which permits 
no additional strain in the spring when the displacement of the end of 
the horizontal member exceeds el. If the end of the horizontal member 
moves upward after exceeding ej_, the ratchet Immediately catches. The 
combined force -displacement relation of each spring system is shown in 
figure 1. This relation may be regarded as a simple stress-strain curve 
that includes a plastic region and the phenomenon of strain reversal. 

If such a column is assumed perfectly straight and Its load Is free to 
change during bending, an infinity of loads can be found at which the 
top of this column can start to assume a deflected position. This range 
of loads is included between the tangent -modulus load for this column, 
and the Euler load. Those loads between the tangent -modulus load and 
the reduced-modulus load require an increase in load during initial 
bending; whereas, those loads between the reduced-modulus load and the 
Euler load require a decrease in load during the bending process. 

The actions of the spring systems at the start of bending can be 
described by locating the instantaneous center of rotation of the rigid 
column. At the tangent -modulus load, the center is at the end of the 
horizontal leg on the side that Is not loading; therefore, no displace- 
ment of that spring system occurs. For loads between the tangent-modulus 
load and the Euler load, the center lies on the horizontal leg between 
the two ends . At the Euler load, the instantaneous center is at the 
opposite side and no displacement of that spring system occurs. It is 
evident, then, that reversal of the spring system on one side always 
occurs for these loads except at the tangent -modulus load when there is 
no change in the displacement of that side. 
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An analysis made for the idealized H-section (appendix C) reveals 
a similar but somewhat more complex behavior of a perfectly straight 
column; that is, the same range of loads exists at which a perfectly 
straight plastic column can start to assume a deflected shape. (See 
fig. 2.) Those loads between the tangent -modulus load and the reduced- 
modulus load are associated with an increase in load during initial 
deflection; whereas, those loads between the reduced-modulus load and 
the Euler load are associated with a decrease in load. The essential 
difference in the column actions for this range of loads is the pattern 
of strain reversal which occurs during bending. At the tangent -modulus 
load, none of the strains over the entire volume of the column reverse 
at the start of bending; however, the strain is Just stationary at the 
center of the column on the convex side. At the reduced-modulus load, 
the reversal is complete over one side of the column as usually presented 
in the double -modulus theory and the strains are stationary along the 
surface that separates the reversed and unreversed regions. At the 
Euler load, all the strains reverse over the entire volume of the column 
except for the center of the column on the concave side, and there the 
strain is stationary. Because the distribution of stiffness over the 
length of the column is the same at these three loads, the instantaneous 
deflected shape is a half sine wave. 


Since the deflected shape of the column at the tangent -modulus and 
Euler loads is known and since the instantaneous center for strain of 
the center cross section of the column is on the convex side at the 
tangent -modulus load and on the compression Bide at the Euler load, the 
initial slope of the curve of load against center deflection can be com- 
puted for the columns. The slopes obtained here are valid for columns of 
constant symmetrical cross section. At the tangent -modulus load. 





and, at the Euler load, 


dP 

d y c 



where b is the distance between the extreme fibers in bending and p 
is the radius of gyration of the cross section. At the reduced-modulus 
load, this slope is zero. 

The analysis of a perfectly straight plastic column leads to a 
range of critical. loadB, but only the smallest, the tangent -modulus load, 
can be accepted as being significant for real columns. No real column 
is perfectly straight; therefore, it is reasonable to define the signifi- 
cant critical load as one based on the behavior of a slightly bent column 
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as the initial lack of straightness vanishes. In order to demonstrate 
the results of such a point of view for defining the critical load, a 
more complete analysis of the spring-supported column was made in which 
a small initial deflection do of the top of the column was included. 

In making such an analysis, it is necessary to keep track of the dis- 
placements of the individual spring systems to be certain that the cor- 
rect force -displacement relations are being used. The possible combina- 
tions for the spring systems are shown in figure 3- 

The results of the analysis of the column with various amounts 
of initial deflection are shown in figure 4. Also shown are the 
limits of the regions in which the various force-deflection relations 
exist. The most significant point to be observed is that, as the initial 
imperfection decreases, the tangent -modulus load is the limiting load 
at which there is a sustained increase in the bending deflection of the 
column. It should also be noted that the reversal of the spring deflec- 
tion always occurs below the tangent -modulus load and occurs just at 
that load as the initial imperfection vanishes. 


THE MAXIMUM LOAD 


The load-deflection curves for the simple column, given in figure k, 
are all approaching the reduced-modulus load at large deflection. This 
is a consequence of the linearity of the spring systems and is to be 
expected if one considers the deflected column with reversal as a new 
elastic column, the stiffness of which is measured by the reduced modulus. 
This behavior is not typical of the actions of real columns, because as 
the load on the column increases, the strains increase and there is a 
continual reduction of the tangent modulus. How rapidly the tangent 
modulus decreases depends on the shape of the stress-strain curve; 
therefore, the maximum strength of a column of a given geometry is 
expected to depend on the stress -strain curve of its material. 

To study the effect of the shape of the stress -strain curve on the 
maximum strength of a column, an analysis was made of the behavior of 
the H-section col umn after it had become critical and started to bend 
at the tangent -modulus load. The analytical stress-strain curve for .the 
material of this column was assumed to be of the form suggested by 
Ramberg and Osgood (reference 7)- This form is summarized in figure 5* 

The stress crp is usually close to the yield stress of the material. 

The most significant parameter for this study is the exponent n. Low 
values of n correspond to gradually curving stress -strain curves, and, 
as n increases, the curvature changes more rapidly at the knee. Values 
of n in the neighborhood of 10 are associated with aluminum alloys; 
whereas, values between 3 and- 5 apply to stainless steels. Magnesium 
and the low-carbon steels have stress-strain curves which correspond to 
values of n of 30 6r greater. 
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The study of strain history of the flanges of a column which starts 
to bend at the tangent -modulus load shows the same general strein history 
regardless of the shape of the stress-strain curve. At the start of 
bending, reversal of stress begins at the center of the convex side of 
the column. As the deflection increases, the region of reversal spreads 
rapidly over the convex side of the column and is complete over the whole 
convex flange at maximum load. After the deflection corresponding to 
the maximum load is exceeded, reversal of strain spreads into the con- 
cave side of the column. 

A summary of the load-deflection results that were obtained for the 
H-section column when the stress-strain curve and the tangent -modulus 
load were systematically varied is given in figure 6. The loads are 
given in terms of a load P]_ that produces an average stress cri in 
the straight column. Two results are significant: First, the smaller 

the value of n, the higher the maximum load in relation to the tangent - 
modulus loadj second, the smaller the value of n, the greater the deflec- 
tion at which the maximum load occurs. The results obtained for the 
maximum load are summarized in two forms: In figure 7> the ratio of the 

difference between the maximum load and the tangent -modulus load to the 
difference between the reduced-modulus load (had the column remained 
straight) and the tangent -modulus load is plotted as a function of the 
tangent -modulus load. The maximum load of columns, critical in the 
plastic stress range, exceeds the tangent -modulus load by a fairly con- 
stant percentage of the difference between Pr^ and P^. In figure 8 
the ratio of Pmax "to Ft is plotted as a function of the tangent -modulus 
load. In the plastic range the percentage increases in maximum load 
over the critical load are roughly constant for a given value of n and 
are larger the smaller the value of n. 

The dashed parts of the curves in figures 7 and 8 have no practical 
significance. Their shape is a consequence of the fact that the Ramberg- 
Osgood stress-strain curves have no proportional limit and are therefore 
nonlinear in what would normally be the elastic range. A more correct 
interpretation would be to consider the ordinates to be unity in the 
dashed regions . 


CONCLUSIONS 


The theoretical study of column behavior in the plastic stress range 
led to the following conclusions: 

1. If the behavior of a perfectly straight column is regarded as 
the limiting behavior of a bent column as its initial imperfection 
vanishes, the tangent -modulus load is the critical load of the column - 
that is, the load at which bending starts. 
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2. The max imum load that an initially straight inelastic column 
can support is larger than the tangent -modulus load. The ratio of the 
max imum load to the tangent -modulus load will he larger for columns 
having stress-strain curves represented by lower values of the 
exponent n in the Ramberg-Osgood representation. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , October 1 6 , 1950 


t 
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APPENDIX A 
SYMBOLS 


ar 

A 

b 


coefficients in sine series 
cross-sectional area 
column thickness 
constants of integration 


do 

d 

e 

e l 

E 


initial lateral deflection of spring- supported column model 
lateral deflection of spring-supported column model 
displacement or strain 

displacement corresponding to force Pi/2 or strain 
corresponding to stress cTj. 

Young's modulus 


®t 


tangent modulus (d a/de) 


Eiji — 


F 

I 

k 




K 


L 


do 

de 


2 


at stress corresponding to tangent -modulus load 
force 

moment of inertia 
proportionality factor (x,= kL) 
spring constants 

parameter used in Ramberg-Osgood representation 
measure of length 


m 

n 

P 


positive integer 

positive exponent in Ramberg-Osgood representation 
variable load on column 
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*1 

r 

u 

x 

X' 

*1,X 2 

y 

y* 

yi^2 

p 

a 

a l 

Subscripts: 

c 

( ^rev 
E 

max 

EM 

T 

L 

R 


load on straight column causing a strain ei 
positive integer used in sine series 

vertical displacement of pivot point of spring-supported 
column 

variable distance along length of col umn 

nondimensional distance along length of column (x = x'L) 

variable dimensions of length measured from particular 
points 

variable deflection 
nondimensional. deflection (y = y’b) 

deflections at xj_ and x 2 
radius of gyration 
stress 

0.7E secant yield stress 
mi&heigdrt of column 

particular value when reversal occurred 

Euler 

maximum 

reduced modulus 
tangent modulus 
left 
right 
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APPENDIX B 

ANALYSIS OF SPRING -SUPPORTED COLUMN 


Two analyses were made of a simple mechanical model of a column. 
First, an infinitesimal-deflection analysis which showed that, if one 
assumed a column to he initially straight, an infinity of equilibrium 
paths could be found by which the column can start to deflect. Each 
path was associated with a definite load. The second analysis was a 
finite -deflection analysis. of the same column and showed that, for real 
columns (which always contain some small initial imperfection), only the 
lowest of these loads, the tangent -modulus load, had any real 
significance . 

If such a column is analyzed by the usual approach, the Euler load 
is 


PE = (k x + kg)^- 


where the spring system stiffness + &2 represents the original 
Blope of the force -displacement diagram. Above the knee of the force- 
displacement diagram, the stiffness of each spring system for increasing 
displacements is kq, and for this column the tangent -modulus load is 


P T 


kib 2 

L 


The double -modulus theory gives for this column a reduced stiffness 
2k 1 (k 1 + k 2 ) 

— r: - and therefore a reduced-modulus load 

2k^_ + k 2 


2ki(ki + k2) b 2 _ ,2 PePt 
P RM = 2^ + 1^ L Pe + P T 

In the development of the analysis it will be convenient to introduce 
these loads. 


Infinitesimal-Deflection Analysis 

Figure 1 shows the details of the model of the column. The vertical 
leg of length L/2 and two horizontal legs each of length b/2 are 
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considered rigid and rigidly Joined. The column is loaded at the top 
of the vertical leg and is constrained at the intersection of the two 
legs so that this point is free to move vertically and rotate but can- 
not move horizontally. At each end of the horizontal member, there is 
set of supporting -springs . In each spring system, the spring with 
stiffness kq is elastic and its grounded end is rigidly supported; 
whereas, the spring with stiffness k 2 has its grounded end attached 


a 


to a ratchet that 


Pi k 2 

slides for a constant force -yr- z 

2 kq 


in the spring. 


+ kg 

On reduction of the force in this spring, the ratchet ceases to slide 
and unloading of the spring occurs with changes in force related to 
changes in displacement according to the stiffness k 2 . The action of 
these two springs of stiffness k]_ and k 2 combine to form a force- 
displacement relationship for each system similar to an idealized stress- 
strain curve with a sharp knee at the force Pq/2 (see fig. 1). 


Let us consider possible changes in the equilibrium position of the 
straight column that has been loaded with some force P which is 
greater than Pq. Static equilibrium of vertical forces and moments 
about the pivot point requires the changes in force in the springs to 
be related to the applied load and to any change in the applied load by 
the following equations: 


AFq + AF R = AP 


AFl - AFr = 2 ^(P + AP) 


> 


(Bl) 


The lateral deflection Ad is assumed to be directed to the left. 
The displacements of the sets of springs can be related by geometry to 
vertical displacement of the column at the pivot point Au and the 
lateral deflection of the top of the column Ad. The displacement of 
the set of springs on the left is 


and on the right. 


Aeq = Aa + Ad £• 


Ae R = Au - Ad £ 


(B2a) 


(B2b) 


In order to complete the analysis of the column, force -displacement 
relations must be introduced for the spring systems. The form of these 
relations depends on three possible displacement patterns: 
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(a) 

eL 

and 

eR 

increasing 

00 

e L 

and 

e R 

decreasing 

(c) 

e L 

increasing and e^ decreasing 


In the following sections, these three displacement patterns are discussed 
and also included in each analysis is the possibility of e^ unfl eR 
not changing. 


Left and right displacements increasing .- The first possible force- 
displacement combination corresponds to the left and right displacements 
increasing or remaining unchanged; that is. 


-AfiL = Au + Ad £ ^ 0 

JLa 

Aeg = Au - Ad jr ^ 0 



(B3) 


For a load P greater than Pi, these changes in displacement multiplied 
By ki give the changes in load as follows: 


AF^ = ^Au + Ad i^ki 
AFr = ^Au - Ad ^kq 


(B4) 


Substitution of equations (B4) in equations (Bl) gives the static equi- 
librium relations that must exist between the change in the column 
deflections and the load: 


2kiAu = AP 


(B5a) 


Ad(pr - (P + AP)] = 0 (B5b) 

kib^ 

where Pqi = — ^ — as noted previously. 

The problem now is to find solutions to equations (B5) which do not 
violate the conditions on the displacements given by the inequali- 
ties (B3). A trivial solution is that Ad is zero and that Au 
increases for a positive AP. Physically, this solution would mean that 
the column remains straight and continues to compress under increasi ng 
load. If, however, in the moment equilibrium equation (B5b), the 
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"bracketed term is zero, it would seem that any Ad is possible and the 
column can bend. During the bending process, however, the load on the 
column. P + AP must always be equal to P^ in order to satisfy moment 
equilibrium} AP must therefore be zero. If AP is zero, Au is 
zero, and only zero Ad can satisfy the inequalities (B3). 


Left and right displacements decreasing .- The second possible 
force-displacement combination corresponds to the left and right dis- 
placements decreasing or remaining unchanged} that is,. 


Then 


Ae T = A 11 + Ad =- ^ 0 


Aeg = Au - Ad =- ^ 0 




£F l = (au + ad + 



Static equilibrium requires that 

2(kq + k 2 )Au = AP 
Ad[pE - (P + AP)] = 0 


(B6) 


(B7) 


(B8a) 

(B8b) 


where Fg = (kq + k2)-^. If AP Is negative, these conditions are 

satisfied for any P If Ad is zero. Physically, this solution 
would mean that the straight column can lengthen for a decrease in 
load. In the moment equilibrium equation (B8b) , it would seem that 
any Ad Is possible if the load on the column P + AP is Pg. 

If, during the bending process, the load P + AP must always be 
Pg, then AP must be zero. If AP is zero, there is no Au, and 

only zero Ad can satisfy the Inequalities (b6) . 


Left displacement increasing, right displacement decreasing . - The 
third force-displacement combination corresponds to the left displace- 
ment increasing and the right displacement decreasing or either 
displacement remaining unchanged} that is. 
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= Au + Ad ^ 0 
Asr = Au - Ad ^ 0 


% 



(B9) 


Then 

AF l = (&i + Ad ^k-j_ 

AFr = ^Au - Ad ^)( k l + ^ 2 ) 

Static equilibrium requires the following relations between the load, 
its change, and the changes in column displacements: 



(2kq + kg)Au - k 2 Ad. ^ = AP 

-k 2 Au + (2ki + kp)Ad — = 2 — (P + AP) 

Lb ) 


(Bll) 


When the previously mentioned relations for Pj, Pg, and Prm are used, 

the displacements obtained from the static equilibrium equations (Bll) 
are 


AuL 



Ad 

b 




r ( B12 ) 




These relations between the change in load and the change in col umn 
deflections are valid only if they are consistent with the original 
assumptions for the force -displacement relations of the spring systems . 
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Substitution of these column deflections in the 
the direction of spring displacement yields the 


Inequalities that define 
inequalities 



(B13) 


Solutions consistent with lateral deflection exist among the equili- 
brium equations (B12) and the Inequalities (B13) when the load on the 
column lies between Prp and Pg. When the load lies between Prp 
and Pgjf, lateral deflection is possible for increases in loadj whereas, 
when the load lies between and Fg, lateral deflection Is possible 

only for decreasing load. If 

*e£*£*bk 

then 

AP 


b 

and if 


then 


P RM = p = P E 


AP 

Pip 

— Ad ~ " 
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Finite -Lateral -Deflection Analysis 
with Small Initial Deflection 


The previous analysis has indicated that an infinity of equilibrium 
paths can be found by which a perfectly straight column can undergo 
lateral deflection. A more realistic approach is to consider an actual 
column as one which admits a small but finite initial imperfection. It 
is of interest, therefore, to analyze the load-deflection history of the 
spring-supported column when such an imperfection exists. 


Consider the same spring- supported column which has as an initial 
imperfection a small deflection do of the tip of the column. Static 
equilibrium requires that 


F L + F R = P 


Fl - Fr = 2P 


do +-d 
b 




(B14) 


Geometry requires the following relation between the spring displace- 
ments and the column deflection: 


e L - e R = 


2db 

L 


(B15) 


The force-displacement relations of the spring systems, as in the 
previous section, depend on the magnitude of the displacements and the 
direction in which they are progressing. During the initial loading, 
the spring systems are elastic. As the load Increases, the spring system 
on the side with the more rapidly increasing displacement becomes 
"plastic" and the column is considered "elastic-plastic." If the initial 
deflection of the column is sufficiently small, the spring system on the 
side with the more slowly increasing displacements can also become 
plastic as the load increases. The region in which this occurs is con- 
sidered "plastic." If, after both sides have become plastic, the 
trailing spring displacement reverses, it does so elastically and the 
column Is considered "plastic with elastic unloading." A summary of 
these possible force -displacement combinations is given in figure 3* 

The column load-deflection relations in these four regions are found as 
follows : 


Elastic range .- In the elastic range the force-displacement rela- 
tions for the springs are 


F L = (^1 + 4 2 )e L 

Fr = (hi + 


(Bl6) 
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Substitution of these relations in the static and geometric relations 
yields the following solution for the lateral deflection of the column: 


_P_ £0 
d _ Ft ~b 
b ~ Pb _ p 

Pip Pip 


(B17) 


where 


PE = 


b g (ki + kg) 
L 


is the Euler load for this column. 


The limit 


of the elastic range is reached when the force in the more rapidly 
loading spring system equals Pi/2. The elastic region in the plot of 
column load against lateral deflection is bounded by the coordinate axes 
and the straight line given by 


d l/Pl _P_ Ft\ . 0 

b 2^ P T Pe j 


Elastic -plastic range .- In the elastic -plastic range a new set of 
force -displacement relations exists for the springs. These relations 
are 


f L - k l e L + fc2 e l 
F r = (k x + k 2 )e R 


(Bl8) 


The lateral deflection of the spring column in this range is defined by 
the equation 


d 

b 



(B19) 


If 



lateral deflection grows in the elastic-plastic range and approaches 
infinity for values of the load approaching as a maximum value the 
reduced-modulus load Prm. Before this max imum load is reached, the 
smaller spring displacement, which has remained elastic, always reverses 
its direction of displacing. For large initial deflections, this 
reversal occurs in the elastic range or at the instant the column enters 
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the elastic -plastic range. If the initial deflections are small enough, 
reversal occurs later. The loci of loads and deflections at which 
reversal occurs in the elastic -plastic region lie on the curve 


d 

b 


+ 




1 _p_ 1 £1 /, 

2 p t + I p T ^ 



if 



the spring displacement on the right, which has always been the smaller 
one, can also become plastic. ,The locus of points at which this change 
occurs is given by the straight line 


d 

b 


P_ 

£t 



= 0 


Plastic range .- The plastic force -displacement relations for the 
springs are 


F L = k l e L + k 2 e l 
F r = k ieR + kg 6 ! 


(B20) 


The lateral deflection in this range is given by the equation 


P d o 


d 

b 


Pt 


_p_ 

Pt 


(B21) 


As the load and lateral deflection increase in the plastic range, the 
spring displacement on the right reverses and, consequently, changes 
the spring force -displacement relation. Reversal of the direction of 
the spring displacement in the plastic range takes place on the curve 
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• Plastic range with elastic unloading . - After reversal has occurred * 
in the plastic range, the force -displacement relations for the springs 
are 


F l = h-|_e L + kgejL 

f R = ( k l + ^2)^ + ^[fl " ( e R)rev] 


(B22) 


The column lateral deflection is defined in this range by the equation 



1 P / P T 

2 f t\ p rm 



P _ lfo 
PRM b 2 \ P RM 





(B23) 


As the load increases, the lateral deflection increases and approaches 
infinity as the column load approaches Prm, the reduced-modulus load. 

Numerical example .- The lateral -deflection analysis of the spring- 
supported column was applied to a particular column. This column is 
defined by choosing a ratio of 

_ 9 

Pi ~ t 

and a ratio of 

£t = 9 

Pp 8 

Such a column has a ratio of 

PRM _ k 
P T 3 

The results of the analysis are given in figure 4 for various values of 
initial deflection of the column. 
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APPENDIX C 

ANALYSIS OF IDEALIZED H -SECTION 


The analysis of the simple spring-supported column, given in 
appendix B, has served the purpose of showing the meaning of the tangent- 
modulus critical load in inelastic column behavior. The model is, how- 
ever, a crude representation of an actual column and has only one pos- 
sible deflected shape. The* analysis fails to reveal another phenomenon 
associated with a real col umn , that is, the growth of the region of the 
column over which strain reversal takes place. A simple model that can 
exhibit the growth of reversal along its length is an idealized H-section 
column. Such a column consists of two flanges of equal area separated 
by a web of negligible thlekness. No shear strain can occur in the web 
so that ordinary beam theory can be used to relate the curvature of the 
lateral deflection of the col umn to the strains in the flanges. 
Furthermore, the state of stress in each flange Is assumed to be one- 
dimensional so that in the plastic range the stresses can be related to 
the strains by the ordinary stress -strain curve for the material. If 
strain reversal occurs In the plastic range, the unloading takes place 
elastically. 


Inf initesimal -Deflection Analysis 


An infinitesimal-deflection analysis of an originally straight 
H-section column reveals the same situation that exists when such an 
analysis Is made for the spring column^ that is, a range of loads exists 
for which a perfectly straight column can start to assume a deflected 
position of equilibrium. 


Derivation of the differential equations .- At any .section along the 
length of the column, static equilibrium requires that 


AFl + AF R = AP 


AFl - AF r = 2 ^(P + AP) 


(Cl) 


Geometry requires that the curvature be related to the strains on either 
side of the column by the equation 


d 2 Ay 

dx 2 


Acl - Acr 
b 


(C2) 
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In order to complete the analysis, the stress-strain relations for 
the flange must be defined. At every section along the length of the 
column there are three possible combinations of directions of straining: 
ep and ep increasing, ex, and ep decreasing, and e^, increasing, 
ep decreasing. Differential equations relating the column deflection 
to the load are derived as follows for these combinations: 

If the strains ep and ep are increasing, the force-strain 
relations for the fl ang e are 


AE t 


^L. - "2 


AEt 

AFp = Aep 


(C3) 


Substituting these relations in the equation of geometry (C2) and 
eliminating AF L and AFp by the use of the static equilibrium equa- 
tions (Cl) yields the following differential equation in which the term 
of higher order has been neglected: 


& 2 A 7 
dx 2 


b 2 AE t 


Ay = 0 


(C4) 


If the strains ep and ep are decreasing, the force-strain 
relations for the flanges are 




= Ae 


AE 

L 


AFp 



(C5) 


Substituting these results In the equation of geometry (C2) anri mBidng 
use of the static equilibrium equations (Cl) yields the following dif- 
ferential equation relating the column deflections to the load: 


d 2 Ay Itf? A 
— o- + - 0 — ^7 
dx 2 b 2 AE 


0 


(C 6) 


If the 
decreasing. 


strains ep are increasing and the strains ep 
the force-strain relations for the flanges are 


are 
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AE+- 

AF l = Ae L -g- 


= ^®R 


7 


AE 

2 


(C7) 


The differential equation relating the increment of column deflection to 
the load and its change is 


a 2 Ay . , 

dx 2 b 2 A\ E t 



£p/_l_ 

Ab^. 


(C8) 


Initial deflection with increasing load . - If the H- sect ion column 
deflects under end load, the sections along its length must he straining 
according to one of the three possibilities previously mentioned. If 
the bending takes place under increasing load, it is reasonable to 
assume that, in the regions near the ends of the column, both strains 
are increasing during bending. At the same time, in a region near the 
center of the column the strains on the left increase while those on the 
right decrease. 

Therefore, differential equation (C 4 ) may be assumed to apply at 
the ends of the column. The solution to this equation is 


Ayi = C-l sin Xq (C 9 ) 

V b AE t 

where Xq is measured from the end of the column. and the cosine term 
has been dropped in order to satisfy the condition of zero deflection 
at the end of the column. At the center of the column, differential 
equation (C8) applies and its solution is 

- c 2 =° s ^(1 + - § b rn| < C10 > 

where X2 is measured from the center of the column and the sine term 
has been dropped because there is no slope at the center of the column. 
Three arbitrary quantities are to be determined: the constants Cq 

and C2 and the increase In load AP. These quantities are determined 
by the conditions that must exist at the junction between the region In 
which all strains increase and the region in which the strains on the 
convex or right side of the column reverse. If the region of increasing 
strains is of length kL, at the junction 
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xj_ = kL 


x 2 


2k 


the following conditions must he satisfied: 

Ayi - Ay 2 = 0 
dAy]_ dAy"2 




dxi 


dX2 


= 0 




= 0 


J 


Substitution of the deflected shapes Ayq and Ay 2 f ’ rom equation 
and equation (CIO), respectively, yields the following homogeneous 
system of equations in C]_, C 2 , and AP: 


sin 


Post 


kPk 2 ! 2 

(b 2 AE t 


|4 Pk 2 L 2 
-AE+ 


-COSI 


J (1 - 2k) 2 L 2 (e + Et) 
2b 2 AEE t 


b_ Et - E 
*2P E t + E 


+ E t (p(l - 2k) 2 L 2 (E + E t ) 
sin\ 


2E 


2b 2 AEE t 


P sin 


[4 Pk 2 L 2 

-AE+ 


b 

2 


Solutions exist for values of the load and regions of reversal (1 
which satisfy the transcendental equation 



(cn) 

(C9) 

Cl 

C 2 = 0 
AP 

(C12) 

- 2k)L 

(C13) 


The principal roots of this equation yield values of load lying between 
the tangent -modulus load and the reduced-modulus load, with corresponding 
regions of reversal, on the convex side ranging from an infinitesimal 
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length at the center of the column to the entire side of the column. 
The shape of the deflected column is a half sine wave at these two 
limiting loads . The rate of change of the load in the col umn to the 
change in deflection at the center of the column is, at the tangent - 
modulus load, 

AP _ 2p T 
A7 c " 

and, at the reduced-modulus load, 


4r c 

Initial deflection with decreasing load .- By a process similar to 
the one used for initial deflection with increasing load, equilibrium 
paths can be .found by which a straight plastic column can deflect under 
decreasing load. Assume that, for a region near the ends of the column, 
the strains on both sides are decreasing; then, differential equa- 
tion (C 6) applies in this region. The solution to this differential 
equation that satisfies the condition of zero deflection at the end of 
the column is 


Ayi 


Co sin 


4P 

b 2 AE 


X 1 


(Cl4) 


In the region at the center of the column, the strains on the convex 
side are assumed to decrease and those on the concave side to increase; 
therefore, differential equation (C8) applies in this region. The 
solution for the deflected shape in this region is 


A72 = c 4 cos 



Xg + 



E t - E 
s + 


(ci5) 


A relation between the arbitrary quantities C3, C4, and AP is found 

from the three conditions that exist at the junction of the region in 
which all strains are decreasing and the region in which only the strains 
on the convex side are decreasing. At the junction, x n = kL and 
1 - 2k x 

x 2 2 ~ — L wJaere 4L is ncfW- the length of the region of decreasing 

strains and the three conditions that must be satisfied are 
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4^1 - Ay 2 = 0 

<3Ay~i dAy 2 
dx-j_ dx 2 

Ae L = 0 


(C16) 


These three conditions lead to the following homogeneous system of 
equations in C 3 , C^, and AP: 


sin' 


jlPk 2 ! 2 

lh 2 AE 


-cos 


P(1 - 2k ) 2 L 2 (E + Et) 
V 2b 2 AEEt 


b_ Et - E 
2P E + Et 


c 3 


COS| 


j kPk 2 L 2 
V h 2 AE 


,F+lt P(1 - 2k ) 2 L 2 (E + Et) 
\— cos> 

' ^t V 2b 2 AEE t 


0 


c 4 


= 0 


P sin 


4Pk 2 L 2 

b 2 AE 


0 


b 

2 


AP 


(C17) 

Consistent solutions exist among the quantities Co, C 4 , and AP, 
which also satisfy the assumed directions of straining if the load and 
the regions of reversal satisfy the equation 


x kPk 2 L 2 ^ |?( 1 

tan \ — ^ tan 


b 2 AE 

and if AP is negative. 


- 2k ) 2 L g (E + Et) _ E + Et 
2b 2 AEE t 1 ^t 


(C18) 


Solutions of these equations exist for loads lying between the 
reduced-modulus load and the Euler load. At these two limiting loads, 
the deflected shapes are sinusoidal. At the reduced-modulus load, strain 
reversal occurs over the whole convex side of the column j whereas, at 
the Euler load, it is complete over both the convex and concave sides. 

At the reduced-modulus load 
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and at the Euler load 


AP 

4/ r 


2P e 


Initial load-deflection relations at tangent-modulus load and 
Euler load for arbitrary cross section .- A general result which can be 
induced from the study of the spring- supported column and the H-section 
column is the initial slope of the load-deflection curve for any 
constant -section column at the tangent -modulus load, at the reduced- 
modulus load, and at the Euler load. At the tangent -modulus load, the 
strain at the center of the convex side of the column is stationary) 
therefore, the instantaneous center of rotation for strains in the cross 
section is at this point. The distribution of strain across the column 
is uniform with a maximum value at the concave side. If the strain on 
the concave side is (AeL ) c , the increase in load would be as follows 

for a symmetrical cross section: 

('Ae L ) c AE T 


AP = 


(C19) 


where A is the area of the cross section. The change in strain Aej, 
is related to the curvature of the column by equation (C2) and is, at 
the center, 


(^ L )c 


•|i (iy sin f ) 


-ay c 


2!_ 

L2 


X= 2 


(C20) 


Elimination of ^A6 l] c between equations (C19) and (C20) yields 


^T * 2 b 
A? = -2" Ay c — 


or 


47 c = 2p 2 

At the reduced-modulus load, the slope is zero. At the Euler load, 
the instantaneous center for strains on the cross section is at the 
center of the column on the concave side; therefore, by a process 
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similar to that employed at the tangent -modulus load., this slope can he 
shown to he 


AP h 

4r c 2p 2 E 


Finite -Deflection Analysis 

An exact finite -deflect ion analysis in the plastic range is com- 
plicated by two factors. First, the deflections may become large enough 
that the more exact definition of curvature of the column is required. 

In the analysis which is developed here, the usual small-deflection 
definition for curvature is used and is justified by the results that 
show the most significant actions to take place for, real columns at small 
displacements. Second, as lateral deflection proceeds., the strains on 
the convex side reverse their direction of straining and elastic 
unloading of the column elements takes place. This second phenomenon, 
referred to as the growth of the reversed region, was not so well under- 
stood and required some investigation. In order to investigate the 
effect of the growth of reversal, a finite-deflection analysis was made 
for an idealized H-section column made of material having a constant 
tangent modulus in the plastic range. Another analysis was made in 
which a more realistic stress-strain relation for the material was used 
to study the effect of the shape of the stress -Btrain relation on the 
strength of columns. 


Analysis of H-sectlon column for constant tangent modulus . - The 
previous analysis of the spring column has shown that a perfectly 
straight column may be assumed to start deflecting at the tangent -modulus 
load. The static relationship between the loads in flanges of the 
H-section and the applied load are given by the following equations: 


f l + F R - p "] 



(C21) 


The geometric relation between the strains in the flanges and the column 
deflection is 


d- 2 ? e L " e R 
dx 2 


b 


(C22) 
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Furthermore , in some part of the column, the strains are assumed to 
continue to increase) then, for constant tangent modulus above Prp, 







(C23) 


where e^ is the uniform strain that exists everywhere in the column 
at the tangent -modulus load. The differential equation which relates 
the deflection to the applied load is obtained by substitution of the 
static and load strain relations into the geometric relation and 1 b as 
follows : 

d 2 y 4 p 

— p + -p — y = 0 (C24) 

dx b AE t 

If at some section of the column the strains on the convex (right) 
side start to decrease, then 


F R 



(C25) 


where the subscript "rev" refers to the load and strain at which 
reversal of strain has taken place. The differential equation relating 
load to deflection along this part of the length can now be written as 
follows : 


+ 2 P ./1 + 

dx 2 b 2 A\ E 



1 P 
b A 




(C26) 


These equations may be put in nondimens ional form by introducing the 
following notations: 




it 2 E<pb 2 A 

Pt * —>aF~ 


P = rt 2 b g A ^t 

RM E + E t 


, = n 2 Eb 2 A 
E 4 l 2 
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In the unreversed region 


a 2 y' + 
d(x ! ) 2 



= 0 


and in the reversed region 


d% 1 
d(x’ ) 2 


+ if- 



ze£ _L 

2 P T 


“(*! 

Irev 

/PT 

P T I 

\ P RM 


(C27) 


(C28) 


The difficulty involved in determining exactly the deflections arises 
from the fact that the growth of the reversed region is not known 
initially but must be determined as deflection proceeds. An approximate 
solution can be made which is based on the method of collocation and 
can be extended to any required degree of accuracy. The solution is 
developed by assuming that the deflections of the column can be expressed 
as the finite trigonometric series of odd sine terms 


m 

y* = ) a r sin(2r - l)itx' (C29) 

Such a series satisfies the boundary conditions and is symmetrical with 
respect to the center of the column. The coefficients of the series are 
obtained from the set of simultaneous algebraic equations obtained by 
substitution of the series into the differential equations (C27) and 
(C28) and satisfying them at m equidistant stations in the half-length 
of the column. Of the two differential equations, the one that applies 
at a given station depends on the direction of straining. It is there- 
fore necessary to keep a check on the direction of straining at each 
station and to use the applicable equation when the direction of 
straining changes. The analysis can be started because the infinitesimal- 
deflection analysis has shown that, at the start of bending at the 
tangent -modulus load, only the strains at the center of the column have 
reversed. As the load increases, the reversed region spreads from the 
center of the column toward the ends . 

The simplest approximate solution is obtained by the use of a 
single, sine term which satisfies only the conditions at the center of 
the column. If 


y’ = y' sin itx' 
c 

substituting this deflection shape in the differential equation (C28) 
and realizing that reversal occurs immediately at the center yields the 
following approximate solution for the deflection of the center of the 
column: 
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The deflection of the center line of the column as P exceeds Pip is 
given in figure 9* Also shown are the results obtained by again 

assuming Ep = '^E with two, three, and five sine terms for the deflected 

shape of the column and by satisfying the equilibrium and stress-strain 
relations at the corresponding number of equidistant points in the half- 
length of the column. All these solutions indicate that the maximum 
load on the column will approach the reduced -modulus load for large 
deflection because of the assumed linearity of the stress-strain curves. 
The solutions further indicate that relatively few terms in the series 
are required for an accurate solution of the load-deflection history 
of the column. 

Analysis of H-section column with a nonlinear stress-strain curve . - 
The finite-deflection analyses of both the spring-supported column and 
the H-section column have indicated that, for a constant tangent modulus 
in the plastic range, the maximum load for a column approaches the 
reduced-modulus value at a large deflection. In order to make a study 
of the effect of a continually decreasing tangent modulus on the strength 
of an H-section column, a finite-deflection study was made by using a 
more realistic stress-strain curve. The Ramberg-OBgood representation of 
the stress-strain curve (reference 7 ) was chosen for the study because 
it fits closely a wide class of materials used in aircraft construction. 
The Ramberg-Osgood representation of the relation of stress and strain 
states that the strain is 


* ■ & + K (if ^ 

where E, K, and n are chosen so that a suitable fit is obtained to 
the stress-strain curve. A recommended procedure for determining the 
three parameters is suggested by Rsmberg and Osgood. The quantity E 
should be the initial elastic modulus, and if 



then Op is the stress at the intercept of the stress-strain curve and 
a line through the origin with a slope 0.7 of the elastic modulus. The 
stress op is approx ima tely the yield stress defined by the 0.2-percent 
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offset. The quantity n is a measure of the sharpness of curvature of 
the stress-strain curve at the yield stress. Smoothly varying stress- 
strain curves that are typical of stainless steel have values of n 
between 3 and 5 . The values of n for the aluminum alloys axe about 10. 
Very large values of n are associated with the magnesium alloys and 
the mild steels. 

The static and geometric relations are the same as those used in 
the analysis of the H-section with constant tangent modulus. Because 
of the form of the Ramberg-Osgood stress-strain relationship, strain in 
terms of force is more convenient. If the strains are increasing, then 



(C32) 


where Pq is the load on a straight column for which the average stress 
is crq, the yield stress. Substituting these strain-force relations 
into the geometric and static relations and making use of the following 
relations 


b L ‘T kl 2 

yields the following differential equation: 


P T “ 


jt 2 Eib 2 A 


I -1 + In& 


Eijt 


7 VPi 


n-l 




d(x’) 2 Pt Wpt\ 
p i T WJ 


sfc y - - ^'> n - ^ - 2 y'> n ] 


The strain-force relations when reversal has taken place are 

3 P 1( ,2F lV ^ 

e L - AE + 7 jffilPq 


e R 


2 Fr 3 Pi 

AE 7 AE 


■(Pr) 


n 


rev 


Pi 


(C33) 


(C3^) 
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These relations give the following differential equation for the 
reversed region: 



(C35) 


After the maximum load in the column is reached, there are sections of 
the column for which the strains on both sides of the column are 
decreasing. The strain-force relations in this region are 


2P l 


S Nrev 

n 

AE 

7 AE 

L ^ 


2P R 

3*1 

PC*)™ 

n 

AE 7 AE 

Pi 

> 


(C36) 


The differential equation which is applicable in this region is' 


d^ ' 
d(x ' )‘ 



ski 


-W 


T 

rev 


Pi 


■M 


-ini 


rev 


Pi 


= 0 (C37) 


The solution of the system of differential equations was carried out by 
the same procedure used in the analysis of the H-section with a constant 
tangent modulus. A solution was assumed in the form of a finite series 
of odd sine terms. Substitution of the series into the appropriate 
differential equations and satisfying them exactly at an equidistant 
number of stations along the length of the" column - yields a system of 
algebraic equations. The system of equations is nonlinear because the 
moduli are a function of the displacements. 

' The simplest approximate solution is obtained by substituting 

y' = y' c sin nx' 

in the differential equation (C35) and satisfying the equation at the 
center of the column. The nondimensional relation between the center-line 
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displacement and the applied load is then obtained from the equation 


(1 + 2y 





28 , 
-y'c- 


PT 

Pi 


3 *| 

7 1 




(C38) 


in which use is made of the fact that reversal occurs at the center of 
the right side of the column at the tangent-modulus load. Solutions of 
equation (C 38 ) for various values of n and for various ratios of the 
tangent -modulus load to the yield-stress load Pq are given in figure 6. 
The same nondimens ional results are obtained in the analysis of the 
spring-supported column. 


In order to study the error due to fixing the shape as a single 
sine curve, solutions were made with three sine curves. These solutions 
were made for n = 2 and 10 and for a ratio of tangent -modulus load to 
yield- stress load of unity. Figure 10 shows a comparison of the results 
of these calculations with those obtained, for a single sine curve. The 
comparison indicates a slight reduction in the maximum load and an 
increase in the deflection at which maximum load occurs. 


The amount of the difference between the tangent-modulus load and 
the reduced-modulus load, which the H-sectlon column could actually sup- 
port were obtained from the analysis summarized in figure 6 and are 
given in nondimens ional form in figure 7 for various values of n. As 
increases, less of the difference is developed for columns critical 
in the plastic range. 

The relation of the maximum load to the tangent -modulus load is 
su mmar ized in figure 8 for various values of n and for columns which 
are critical at various stress levels. 

In figures 7 and 8, the dashed parts of the curve were not calcu- 
lated and are estimates of what would be obtained from the analysis 
using the Ramberg-Osgood form of the stress-strain relation. These 
dashed parts of the curve are trivial and. misleading and are caused by 
the fact that the analytical form of the stress -strain curve has no 
truly linear elastic region. In the elastic region, no difference- 
should exist between the ‘maximum load, the tangent -mo dulus load, 
and the reduced-modulus load (large deflection effects excluded). 
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o Left spring system 
□ Right spring system 



(a) Elastic. (b) Elastic- plastic. 



(c) Plastic. (d) Plastic with 

elastic unloading. 

Figure 3*- Four force -displacement combinations for the spring systems. 
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Figure 4(a).- Load-deflection curves for spring-supported column model with 

initial deflection. 
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Figure 6.- 
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Figure Load-deflection curves for idealized H-aection column with 
constant tangent modulus in the plastic range. 






